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ABSTRACT 

We present the first investigation of Th abundances in Solar twins and analogues to understand the possible 
range of this radioactive element and its effect on rocky planet interior dynamics and potential habitability. The 
abundances of the radioactiv e elements Th and U are key components of a planet’s energy budget, making up 
30% t o 50% of the Earth’s (lKorenaga| 2008 ; Allegreetak |2001b lSchubert et al.ll98(llLvubetskava & Korenagal 
120071 IThe KamLAND Collaborationl201 LlHuang et al. 112013h . Radiogenic heat drives interior mantle convec¬ 
tion and surface plate tectonics, which sustains a deep carbon and water cycle and thereby aides in creating 
Earth’s habitable surface. Unlike other heat sources that are dependent on the planet’s specific formation his¬ 
tory, the radiogenic heat budget is directly related to the mantle concentration of these nuclides. As a refractory 
element, the stellar abundance of Th is faithfully reflected in the terrestrial planet’s concentration. We find 
that log e T h varies from 59% to 251% that of Solar, suggesting extrasolar planetary systems may possess a 
greater energy budget with which to support surface to interior dynamics and thus increase their likelihood to 
be habitable compared to our Solar System. 

Subject headings: planets and satellites: terrestrial planets, planets and satellites: tectonics, planets and satel¬ 
lites: interiors, planets and satellites: physical evolution 


1. INTRODUCTION 

The Earth is a dynamic planet, in which the “heat engine” 
that powers the present-day dynamics that include plate tec¬ 
tonics, the geodynamo in the core, and the ongoing recycling 
and sustentation of the Earth’s atmosphere. Due to melting 
during the planetary formation process, the Earth has dif¬ 
ferentiated into two main compositional parts: the metallic 
core and the silicate-dominated mantle and crust. The metal¬ 
lic core contains iron, nickel, and those siderophile elements 
that partition into the metallic phase. The bulk silicate earth 
(BSE) contains the balance of the Earth’s composition, in¬ 
cluding most of the oxide-forming elements. This includes K, 
Th, and U, whose radioactive isotopes account for 20.1 TW 
(43% . iH uang et a fll2013l) of the Earth’s 47 TW surface heat 
flow dDavies & Daviesll2010l) . The remaining 26.9 TW are a 
result of the initial differentiation into mantle and core, sec¬ 
ular cooling of the core and mantle and the crystallization of 
the inner core. Geodynamo models show that the heat flow at 
the core mantle boundary from the core has remained roughly 
constant at 9 TW (iLabrosse et al.l[l997h . This is ^20% of the 
Earth’s total energy budget today and less than 10% in the 
Archean, including the effects of potassium. 

Within the BSE, however, K, Th, and U are not distributed 
homogeneously. These lithophile elements partition into the 
magma generated at subduction zones, building continental 
crust that accounts of 40% of the Earth’s surface, a direct con¬ 
sequence of plate tectonics. This process leads to heat produc¬ 
tion of 7.6 TW within the continental crust, which serves as 
an insulating layer and thickens the surface boundary layer. It 
is the remaining 12.9 TW of the radiogenic heat budget that 
drives mantle convection. 
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Of the two main planetary heat sources (radiogenic and sec¬ 
ular cooling), the power produced by radionuclides is entirely 
compositionally dependent with the radiogenic heat budget 
inherited upon formation. Secular cooling, in contrast, is de¬ 
termined by the specific formation and segregation history. 
Our goal here is to explore the degree to which the radio¬ 
genic heat budget can vary between terrestrial exoplanets, and 
consequences of formation history will be left to subsequent 
work. 

U and Th are both refractory elements with 50% conden¬ 
sation temperatures of 1604 and 1647 K, respectiv ely, in 
protoplanetary disks of Solar composition (!Loddersll2003tl . 
The measured absolute abundances of other refractory, rock¬ 
building elements in Solar twins and analogues, stars of sim¬ 
ilar metallicity, mass and surface temperature to that of the 
Sun, varies between stars. Furthermore, the Sun is moderately 
depleted in many of the highly refrac tory e lements (Tc on d > 
1500 K) of Z < 56 ( Melendez et a l.l BOP^ . If these differ¬ 
ences are due to intrinsic variations in primordial abundances 
in planetary systems, rather than to the process of planet 
formation, then this suggests that the stellar abundances of 
U and Th may also vary within the Galaxy. In both older 
metal-poor stars and younger dwarfs/subgiants, wide varia¬ 
tion in the absolute abundance of Th has been found, includ¬ 
ing other heavy r-process elements , to track Galactic chem¬ 
ical evolution dMorel et all 1 19921 [Me W illiam et all 119951: 
Johnson & BolteTfeooU iHonda et alJ l2004t Idel Peloso et all 

2005allbh . It is not currently known, however, whether these 
variations in Th or U are present in younger, metal-rich, plan¬ 
etary host stars. Any variation in these abundances will trans¬ 
late into similar increase or decrease in both the crust and 
mantle of any terrestrial planets orbiting these stars thus hav¬ 
ing implications for these heat budget of these planets and 
potential for dynamics. 

2. SAMPLE AND METHODS 
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Table 1 

Input stellar model parameters and measured Fe, Th and Si abundances. References for stellar input model or line parameters are shown as superscripts. 





Fe 4 

Th 5 

# of 

Si 6 


Star 

Age (Ga) 2 T eff (K) log(g) (cgs) 

[Fe/H] 

loge 

loge 

spectra 

(1) (2) (3) (4) 

Average 


with observed planets 


Sun 1 

4.5 

- 

- 

- 

7.47 

0.09 

1 

- 

- 

- 

- 

- 

Sun (This Study) 

4.5 

5777 

4.44 

0.00 

7.47 

0.09 

1 

7.58 

7.60 

7.55 

7.55 

7.57 ±0.02 

HD98649 2 

2.1 

5760P25 

4.51±0.05 

-0.01i0.02 

7.36 

0.31 

1 

7.54 

7.59 

7.54 

7.54 

7.565 ±0.03 

HD1461 2 

4.8 

5756P44 

4.37±0.05 

0.189i0.015 

7.60i0.017 

0.33P0.05 

9 

7.74 

7.79 

7.77 

7.74 

7.76±0.02 

HD16417 2 

7.0 

5812±34 

4.09P0.05 

0.094i0.004 

7.53i0.03 

0.31 ±0.05 

7 

7.64 

7.67 

7.64 

7.64 

7.65 ±0.01 

HD102117 2 

5.6 

5690±22 

4.30±0.04 

0.304i0.003 

7.76i0.02 

0.49±0.04 

8 

7.85 

7.91 

7.95 

7.90 

7.90±0.04 

HD141937 2 

1.3 

5900P19 

4.45±0.04 

0.125i0.003 

7.53i0.02 

0.44±0.02 

3 

7.63 

7.68 

7.73 

7.68 

7.68 ±0.04 

HD 160691 2 

4.6 

5809±22 

4.28±0.04 

0.298i0.003 

7.76i0.02 

0.44±0.09 

9 

7.80 

7.85 

7.95 

7.85 

7.86±0.05 

HD10700 3 

- 

5312±137 

4.59±0.13 

-0.43i0.15 

6.91i0.05 

-0.14±0.05 

4 

7.37 

7.37 

7.22 

7.22 

7.30±0.08 


without observed planets 


HD115169 2 

1.7 

5815±22 

4.52±0.05 

-0.01±0.02 

7.41 

0.36 

1 

7.59 

7.59 

7.54 

7.59 

7.57±0.03 

HD146233 2 

3.1 

5822±9 

4.45±0.02 

0.051 ±0.002 

7.48±0.01 

0.42±0.04 

4 

7.45 

7.60 

7.60 

7.60 

7.56±0.07 

HD45346 2 

3.8 

5745±25 

4.47±0.05 

-0.01±0.02 

7.46 

0.16 

1 

7.54 

7.59 

7.54 

7.54 

7.55±0.02 

HD59711 2 

2.7 

5740±15 

4.50±0.05 

-0.092±0.02 

7.35±0.05 

0.20±0.02 

3 

7.56 

7.46 

7.46 

7.51 

7.50±0.04 

HD71334 2 

4.7 

5724±15 

4.46±0.03 

-0.044±0.02 

7.34±0.04 

0.24±0.03 

5 

7.51 

7.56 

7.56 

7.58 

7.55±0.03 

HD78534 2 

4.6 

5800±25 

4.41±0.05 

0.07±0.02 

7.49 

0.19 

1 

7.57 

7.62 

7.62 

7.62 

7.61±0.02 

HD78660 2 

1.5 

5782±29 

4.52±0.04 

0.033±0.02 

7.45±0.01 

0.26±0.02 

5 

7.63 

7.56 

7.58 

7.58 

7.61±0.03 


References. — 1. Lodders < 2003) 2. Baumann et al. < 20 10) 3. Maldonado et al. < 2012) 4. Sneden et al. (19961) 5. iKupka et all(200CI) 6. Ida Silva et alj(20121) 

Note. — Fe: A = 4019.053 A , log#/ = -2.91, x = 2.609 eV. Th: A = 4019.129 A , loggf = -0.27, x = 0.0 eV. Si (1): A = 5517.533 A , loggf = -2.42, x = 5.08 
eV. Si (2): A = 5621.607 A , log gf = -2.61, x = 5.08 eV. Si (3): A = 5665.563 A , log gf = -1.97, x = 4.92 eV. Si (4): A = 5684.484 A , log gf = -1.65, x = 4.854 eV. 



Figure 1. The spectrum of HD146233 around the Th I line at 4019.129 A 
(dots). The synthetic spectra (solid lines) were computed adopting the stellar 
parameters listed in the text. Synthetic spectra were calculated using log exh 
= 0.17, 0.25, 0.35 and 0.45. The best fit for this spectrum (logexh = 0.35) is 
shown in bold. 

Our sample of fourteen stars (seven with and seven 
without observed planets) was selected from the publicly 
available HARPS archive ( Mayor et al . 12003} ). Obser¬ 
vations were made using the 3.2 m La Silla telescope 
(R = 115,000). In order to minimize systematic uncer¬ 
tainty from input synthetic stellar models^ only known 
Solar twins and analogues from iBauma nn et all (120101) : 
iMaldonado et alJ (l2012l) : lRamfrez et al.l (120091) were selected 
(Table []}. One-dimensional, hydrostatic, plane-parallel stel¬ 
lar model atmospheres were obt ained from the MARCS 
catalogue dGustafsson et al.l 120081) . with raw spectra nor¬ 
malized using the spectrum analysis software SPECTRE 
(Fitzpatrick & Sneden 1987) wit h inpu t stel lar model param- 
eters taken from lBaumann et all (120101) and iMaldonado et all 


(120121) . Th abundances were measured with the line anal¬ 
ysis software MOOG (!Snedenlfl973l) using the blended Fe- 
Th line at 4019.053-4019.129 A adopting the line list of 
ISneden et al.idl996i) supplemented with lines from the VALD 
catalogue (IKupka et all 120001 Table [TJ. We first fit the Fe 
line at 4019.053 A, followed by varying Th until a best fit 
was found (Figure [4} as determined by a minimized \ 2 test. 
Our Fe abundances (TableQ] FigureQ} are in good agreement 
with previous measurements for the same Solar twins and ana¬ 
logues, suggesting that our measurements are accurate with 
respect to Fe. In order to determine if any error is introduced 
by uncertainties in the stellar models, each abundance was 
remeasured adopting values 100 K, 0.05 (cgs) and 0.05 dex 
above and below the published value from the literature for 
T e ff, log(g) and [Fe/H] respectively (TableQ}. These values 
are at or above the published errors for these measurements. 

There are several possible blends in our line list in addition 
to the well-separated Fe I line within this spectral region in¬ 
cluding contributions from Ti, V, Cr, Mn, Co, Ni, Ce, Nd, and 
Sm (Figure|2}. As noted in Section 3, this blending is of par¬ 
ticular importance as it may account for differences between 
our measured Th abundances and those from earlier work 
on similar stars, particularly that of ldel Peloso et al] (12005al) . 
Of these elements, we find that only varying Co made any 
significant difference in our measured Th abundance, while 
changing the abundance of any these other element affecting 
only the fit to the Fe line or required variations well beyond 
the me asured abundances of these elements (del Pelo so et all 
12005al) before the Th abundance was affected. Furthermore, 
two st ars as well as the Sun ov erlap between our sample and 
that of ldel Peloso et alj (l2005ai) . with our line lists being iden¬ 
tical. Adopting their reported Co, Fe and Th values for these 
stars, we are either unable to fit the measured spectra or no no¬ 
ticeable difference was found between their synthetic spectra 
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and ours (Figure 0. Furthermore, Idel Peloso et al.l (12005ah 
fits ~ 8 Co lines for these stars, making their measured Co 
abundances independent of the Th value. We therefore con¬ 
sider any contamination to be minimal and simply adopt the 
scaled-solar MOOG abundances for these elements when de¬ 
termining Th abunda nce. _ 

Not considered by Idel Peloso et ail d2005al) . however, is a 
CH line blended at 4019.170 A. We find that variation in the 
l2 C/ 13 C ratio was found to have a negligible effect on the 
spectral model for this line, and we therefore set this rati o to 
the Solar value of 89.4 for all measurements dAsplund et al.l 
120091) . [C/Fe] for each star was assumed to scale from So¬ 
lar with metallicity as calculated by MOOG and provided a 
good fit to each of the l2 C lines observed. iMelendez et al.l 
d2009t) found that [C/Fe] does not scale with Solar metallicity 
for Solar twins and analogues, however, we find no difference 
in depth the I3 C line when the input, MOOG C abundance is 
scaled by their measured variation of between 0 and 0.06 dex 
from Solar. 

As we discuss in Section 3, Si is the relevant rock-building 
element in terrestrial mantles. To measure Si abundance 
then, we adopt the line parameters for 4 Si lines listed in 
Ida Silva et alj d2012i) . The values for Si abundance derived 
from these individual lines for our sample of 14 stars are listed 
in Table Q] Furthermore, we find only minor variations in de¬ 
rived Si abundance between these individual lines for each 
star in our sample. Again, each abundance was remeasured 
using models of 100 K, 0.05 (cgs) and 0.05 dex above and be¬ 
low the published value from the literature for T e ff, log(g ) and 
[Fe/H] respectively, and no noticeable difference in Si abun¬ 
dance was found. 

3. RESULTS 

For this sample, we find logexh between -0.14±0.04 and 
0.49±0.05fl Our measured Solar log en, of 0.09, in agree¬ 
ment with the accepted literature value dLoddersI 12003 ). We 
find that the Sun is depleted in Th/H relative to 13 stars in our 
sample, with one star more depleted than the Solar abundance 
(Figure0. This corresponds to a variation between 59 (+7/-7) 
to 251 (+24/-22)% compared to the Solar abundance. Our to- 
tal mea sured variation is roughly the same as Idel Peloso et al.1 
( 2005ah . however, in general our measured Th abundances are 
larger. While this enrichment is unlikely caused by our not 
accounting for potential contamination from other elements 
(Figure 0, it may be due to our two datasets varying in both 
resolution (R = 48,000 vs. R = 115,000) and differences in 
continuum normalization. Our measured Solar Th and Fe 
abundances, howeve r, are in agreement with the meteoritic 
value (Table 0, with del Peloso et al. (2005a) measur ing de¬ 
pleted Th and enriched Fe abundances relative to [Codders 


Differences in stellar age within our sample are another pos¬ 
sible source of this variation. We find, however, only minor 
differences in Th abundance when these differences are taken 
into account (Figure0. Furthermore, the half life of 232 Th is 
14 Ga or roughly the age of the Universe, no significant varia¬ 
tion on account of varying ages within our sample is expected. 
Assuming a star’s Th abundance reflects a primordial, homo¬ 
geneous abundance present at the birth of every star in our 
Galaxy, our measured 251% variation cannot be explained. 
The variation in Th abundance we observe here is likely due 

1 log^=log 10 (^)+12. 



Figure 2. The region between 4018.95-4019.25 A including the dominant 
spectral lines within this region as taken from VALD. The Solar spectrum is 
included for reference. 


then to individual stars being seeded with differing amounts 
of r-processed material with these abundances being a lower 
limit on the initial variation. 

A higher Th/H in a star’s atmosphere does not directly in¬ 
dicate a higher Th concentration in a terrestrial planet’s man¬ 
tle. Within the Earth, the majority of Th is present in the 
silicate mantle ([McDonough! 1200 3). Si is also a refractory, 
lithophile element in systems of Solar and near-Solar com¬ 
positions dLoddersI120031) . If we adopt then stellar Si abun¬ 
dance as a proxy for potential terrestrial planet mantle mass, 
normalizing to Th/Si allows us to compare Th concentration 
across planets of different silicate mantle thickness. A planet 
forming from materials with a high Si abundance will, to first 
order, result in planets with a silicate mantle accounting for 
a greater fraction of its radius. With this normalization then, 
absolute Th abundance is a function of planetary mantle mass 
only. The ability for a planet to convect is strongly dependent 
on the thickness of the mantle layer. Small planets forming 
with high Th/Si then may be more apt to convect due to this 
increase in internal heat rather than transport this heat via con¬ 
duction. The Earth’s BSE has an effective logf'l h/Si)~-7.34 
(lMcDonoughll2003l) . a 28% enrichmen t of Th relative to the 
Sun (log(Th/Si)=-7.45, iLodderslfeOOSl) . This may be a con¬ 
sequence of small amounts of Si incorporated into the Earth’s 
core, constrained via joint mineralogical and seismic models 
to be no more than 4.5 wt% (best fit 1.9 wt%). Accounting for 
Si in the core, the Earth’s log(Th/Si) is greater than or equal 
to -7.38 (a 10% decrease). This decrease is small compared 
to our observed 251 (+24/-22)% variation in Th abundances 
and thus normalizing to Si provides us with a good proxy for 
planetary mantle Th concentration. 

The abundance of Si for majority of our sample falls be¬ 
tween loge# of 7.5 and 7.9 including that of the Sun (Figure 
0. Within this range, however, Th/Si varies by a factor of 
2.6, suggesting that the observed variation in Th indeed corre¬ 
lates with an increased concentration of Th within the silicate 
species. We further find a positive correlation of Th with Si 
(c/Th/t/Si = 0.77). We do not observe any variation of Th/Si 
with T e ff and therefore do not suspect any significant selec¬ 
tion effects in our sample. We see this enrichment relative to 
Solar in stars both with and without observed planets, sug¬ 
gesting that although the planet-hosting stars contain known 
gas giants, these enrichments are not indicative of the infall of 
refractory material back onto the parent star. The variance in 
our observations suggest that any terrestrial planets present in 
these systems may be more or less dynamic due to variations 
in their energy budget from radiogenic heat as compared to 
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4019.0 4019.1 4019.2 

Wavelength (A) 


Figure 3 . Syn thetic spectra for the three stars measured both by this study 
an d|del Peloso et all i2005a ). Solid lines represent our synthetic spectra (Ta¬ 
ble!]}. dashed lines are syntheses adopting their reported Fe, Co and Th abun¬ 
dances and dotted is our reported Fe and Th abundances but adopting their 
Co abundance. Squares are the corresponding spectra. 

the Solar System. 


a planet’s dynamic state due to this change in a planet’s 
internal energy budget. This model relies on the sys¬ 
tem’s Rayleigh number as a measure of convective vigor. 
While a parameterized convection model is not as accu¬ 
rate as a full, 3D convection model, it allows us to ex¬ 
plore the general dynamical effects of incorporating a vary¬ 
ing proportion of heat-producing radionuclides into a plan¬ 
etary mantle and to guide further investigation using more 
exact models. There is no mechanism to induce or sus¬ 
tain plate tectonics without interior convection and although 
plate tectonic-like regimes are a more complex question 
of fault st rength, surface gravity, and the presence of liq¬ 
uid water (l\^j encia_&O!Connglj|2()0^;|\^lenciaeta]J|2007|; 
Crowlev et al 1201 HI van Heck & Tackle\l201 iHTacklev et alJ 

20131) , the Rayleigh number provides, to first order, a measure 
of whether interior mantle convection, and thus the possibil¬ 
ity of tectonics will occur in these planets and be potentially 
habitable. 

The Rayleigh number (Ra) of a system as a function of time, 
f, is: 

Ml|= ^° Ar|l|i> ’ c ' (i) 

kvo(T) 

where g is the acceleration due to gravity, po is the average 
density of the mantle, a is the thermal expansivity of the ma¬ 
terials present, A T(t) is the change in temperature across the 
convecting layer, D is the thickness of the convecting layer, C p 
is the specific heat of the system, k is the thermal conductivity 
and 770 is an average viscosity. Above the critical Rayleigh 
number of 10 3 dSchubertill979i : ISchubert et al.ll 1980b . mantle 
interior heat will be more efficiently transported via convec¬ 
tion, with conduction dominating below the critical value. 

A parameterized convection model relates the efficiency 
of heat extraction from the planet to the vigor of convec¬ 
tion, which, in turn, is dependent upon the internal tempera¬ 
ture. The temperature change in a planetary body is quantified 
through conservation of energy, such that: 

MmC P ^- = H (t)-Q(t) (2) 

where H(t) is the amount of heat generated within the mantle, 
and Q(t ) is the heat extracted from the outer surface and M m 
is the mass of the mantle. Heat extraction at the surface is 
limited by the thickness of the surface boundary layer, 6. We 
adopt the relationship: 


3.1. Thermal Effects ofTh budget on Planetary Dynamics 

The nucleosynthetic origins of T h and U are observation- 
ally and theoretic ally correlated dGorielv & Arnouldl 120011 
Frebel et al. 2007) as well as direc tly measured in meteorites 
and crustal minerals i n the Earth (lLoddersll2003l llJrevll 1 952k 
iRudnick & Gaoll2003l) and both must be considered in the to¬ 
tal heat budget of a planetary mantle. If we adopt the aver¬ 
age measured present-day BSE value of Th/U = 4 and scale 
the current BSE Th abundance from geoneutrino studies of 79 
ppb dHuang et alJ2013l) by the maximum Th/Ths un in our data 
set, we find a planet’s energy radioactive heat budget will be 
on average ~ 26 TW greater than the Earth over 12 Ga after 
planetary formation (Figure^. 

3.2. Thermal Model 


6 = D 



(3) 


where /3 is a scaling factor that relates the vigor of convection 
to the surface heat flux, Q(t). The value of (3 for the Earth is 
a matter of debate, but is general considered to be between 0 
and 0.33 and we adopt a median value of 0.2 for this model. 
Equation ([3]) may be rewritten in terms of the surface heat flow 
in Watts as Q(t ) = 47rr 2 * kAT/5: 


Q(t ) = 47rr 2 


kAT(t) 

D 



where r is the radius of the planet. 
Combining equations fll-S yields: 


(4) 


We ado pt a single layer parameterized convection model 
(iMcNamara & van Kekenl 120001) to determine the effects on 


dT 
p dt 


M m C p '-j- =//(7) —47rr 2 


kAT {t) 
D 


/ gP laAT(t)D 3 C p \ 
V kr]o(T)Ra cri , J 


(5) 
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Figure 4. Stars with observed planets are shown as closed squares and without as unfilled squares. The Sun is shown for reference as a triangle. Left: Abundance 
of Th as a function of Fe. The solid line represents the best fit to those stars belonging to the high Fe. high Th group and the dashed represents a one-to-one 
relationship. Center: Abundance of Th as a function of average Si abundance adopting the same legend as in Th vs Fe. Right: Th/Si as a function of average Si 
abundance. For each plot stars with observed planets are shown as closed squares and without as unfilled squares. The Sun is shown for reference as a triangle. 



Figure 5. Same as Figure |4}T with Th abundances corrected to t = 0 Ga 
(circles). Stars without planets are shown as open symbols, stars without as 
closed and the Sun as triangles. 

This captures the basic physics of the feedback loop be¬ 
tween temperature, internal viscosity, and surface heat loss, 
by which the increased vigor of convection increases sur¬ 
face heat loss, thereby reducing internal temperatures, in¬ 
creasing the viscosity and thus limiting internal convection 
(iChristensenli 19851) . This model assumes that any heat trans¬ 
ported to the surface is removed from the system with no in¬ 
sulating thermal boundary layer at the top of the convection 
zone (i.e. continents). While not strictly valid in the presence 
of an insulating atmosphere, this assumption allows us to see 
the first-order effects of an increase in overall thermal con¬ 
ductivity and viscosity profile of a planet with increased Th 
and U abundance. It also does not take into account any feed¬ 
back from tectonic recycling of cooler surface plates into the 
hot mantle, and this model is therefore meant as a first order 
approximation of the thermal history of a planet to determine 
the effects of composition on bulk thermal transport. 

We adopt Earth-like values for C p , k and a (Tabled- In the 
high temperature limit, C p is a constant, and a and k will vary 
by less than an order of magnitude due to bulk mineralogical 
differences and temperature effects, thus expressing very lit¬ 
tle variability in equation 0 In calculating A T(t), we adopt 
A T(t) = T(t)-T s , where T s is the temperature at the top of the 
convecting layer. In our model we hold T s , at a constant tem¬ 


H (TW) Viscosity (Pa s) 



Time After Formation (Ga) Time After Formation (Ga) 

Figure 6. Power produced from radiogenic sources, viscosity, temperature 
and Rayleigh number for a 1 M®, 1 R.|, planet of Earth-like composition 
with 198.3 (red), 79 (black) and 46.6 (blue) ppb Th up to t = 12 Ga. Models 
were run with an initial temperature of 1500 K (dashed) and 3500 K (solid) 
with /3 = 0.2. The approximate critical Rayleigh number (10 3 ) is shown for 
reference as a dotted horizontal line. 

perature of 273 K. 

The thermal evolution model is dependent upon an initial 
average temperature of the planet’s interior (To). This reflects 
the formation history of the planet and is primarily depen¬ 
dent upon two sources of heat: gravitational energy converted 
to heat via both planetary accretion and the segregation of 
a central core. Therefore we assert that varying To serves 
as a proxy for the amount of latent primordial heat present 
via these history-dependent mechanisms, with high 7o corre¬ 
sponding to either planets that formed quickly (and thus re¬ 
tained much of their accretionary heat) or segregated a large 
core relative to their overall volume. 

3.2.1. Mantle Heat Generation 

Heat production, 77(f), arises from the sum of the energy 
produced via the radioactive decay of 235 U, 238 U, 232 Th and 
40 K. We adopt Earth-like mantle heat production as our ba¬ 
sis for scaling Th abundance with present-day values of Th/U 
= 4 and K/U = 10 4 and a BSE value of 232 Th = 79 ppb of 
[Huang et all d2QT3h with 238 U/U = 0.9928, 40 K/K = 1.19* 10“ 4 
(Turcottej 1980f) . These values, as well as the respective ratios, 
are scaled back in time to determine an 77(f = 0) of 51.6 TW 
consistent with the back extrapolation in time from the present 
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Table 2 

Rayleigh number parameters adopted in our model 


Parameter 

Value 

Reference 

r (km) 

6371 

- 

D (km) 

2873 

i 

g (m s“ 2 ) 

9.81 

i 

pa (kg nr 3 ) 

4400 

i 

M m (kg) 

4.04* 10 24 

i 

C p (JK- 1 kg- 1 ) 

1250 

i 

a(K->) 

3*10- 5 

i 

k(Wnr‘ K-‘) 

5.6 

i 

A(s -1 ) 

2.67*10 17 

scaled to diffusivities 

m 

2.5 

assumed same as olivine 

d (m) 

0.005 

equal to xenolith grains 

b (m) 

6*10- 10 

scaled to olivine lattice 

Ai(GPa) 

184 

2 

770 (Pa s) 

3.13*10 9 

calculated 

E* (kJmoL 1 ) 

502.54 

3 

V* (cc moL 1 ) 

2.0 

4 


Referenc es. — 1. [Mc Namara & van Keken i 200fl) 2. 


Keppleretal. 1 2008) 3. Dobson et al. 1 2008) 4. Holzap fel et alj 

goof) 

value of ll(4.5Ga) = 16.4 TW. While a significant fraction of 
Th is concentrated in Earth’s continental crust, this fractiona¬ 
tion is a direct consequence of melting associated with plate 
tectonics and mantle convection modeled here. Therefore, we 
choose to model initially the influence of variation in this pri¬ 
mordial Th abundance on planetary mantle dynamics rather 
than the more complex question of sequestering some of these 
radionuclides in the crust due to plate tectonics. 

3.2.2. Viscosity 

Power-law creep is the dominant creep mechanism in 
Earth’s mantle, which is dependent upon the convective 
driving stress. For simplicity sake, we consider a diffu¬ 
sion creep model (e.g. iKarato & Wu 1993 ) to capture 
the broad temperature and material dependence on viscos¬ 
ity. We assume that perovskite viscosity scales relative to 
dry olivine diffusion creep using the equation for strain rate, 
e = A(cr/n)(b/d)'"exp[-(E*+PV*/RT)], where A is the pre- 
exoponential factor, a is the shear stress, /i is the shear mod¬ 
ulus, b is the length of the Burgers vector, d is the grain size, 
m is the grain size exponent, E* is the activation energy, V* is 
the activation volume, R is the gas constant, P is the pressure 
and T is the temperature. The effective viscosity is then: 

ri(T) = ryo exp [(E*+PV*) /RT} (6) 

where 770 is the viscosity coefficient for a fixed driving stress. 
Viscosity model parameters are as in Table [2] All viscosity 
calculations were calculated at a constant pressure of 86 GPa, 
which roughly corresponds to the pressure at one half of the 
Earth mantle’s radius. 

3.3. Results of Thermal Model 

Solving Equation[5]using a 4th order Runge-Kutta method, 
we find that for planets forming with a high initial tempera¬ 
ture (3500 K), the Rayleigh number remains above the criti¬ 
cal value for convection regardless of Th concentration over 
the entire time period considered here. For planets forming 
with a low initial temperature (1500 K), characteristic of a 
slow coalescence from planetesimals, however, a factor of 
2.51 increase in Th abundance reduces the time at which the 


planet’s Rayleigh number rises above the critical value by ~3 
Ga. When the Th abundance is decreased by a factor of 0.59 
this crossover time increases by 7 Ga compared to an Earth¬ 
like case. These results, however, do not consider the addi¬ 
tional potential variability due to radioactive 40 K, a moder¬ 
ately volatile element with a 50% condensation temperature 
= 1006 K dLoddersl2003l) . While an increase in a planets 40 K 
abundance will further increase the chances of convection in 
a planet, stellar abundance measurements of K will not as di¬ 
rectly inform its contribution to the potential heat budgets and 
habitability of terrestrial planets as compared to the refractory 
elements U and Th. 

4. CONCLUSION 

The measurements presented here are the first measure¬ 
ment of a planetary system’s potential energy budget and po¬ 
tential for planetary interior dynamics. Mantle convection, 
as the underlying mechanism driving plate tectonics on the 
Earth, recycles and regenerates crus t al, oceanic and atmo¬ 
spheric material (ICrowlev et alJl201lt lKorenagall201 ll) . This 
surface-to-interior process regulates the global carbon and 
water cycles, which in turn, aid in creating a habitable sur¬ 
face (iSleep & Zahnlell2001 ). We find that some stellar sys¬ 
tems possess a greater Th abundance compared to Solar and 
due to the refractory nature of Th, this increase will be re¬ 
flected in the mantles and crusts of any resulting terrestrial 
planets. Planetary systems with higher Th/Si ratios at t = 0 
Ga than Solar suggest that these systems possess larger en¬ 
ergy budgets. This supports interior dynamics, and increases 
the likelihood for carbon and water cycling between the sur¬ 
face crust and planetary interior, thus broadening the range of 
planets which may support habitable surfaces. 

This work is supported by NSF CAREER grant EAR- 
60023026 to WRP. 
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